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Abstract 

We propose a discrete traffic flow model with dis- 
crete time. Continuum limit of this model is equiva- 
lent to the optimal velocity model. It has also an ul- 
tradiscrete limit and a piecewise- linear type of traf- 
fic flow model is obtained. Both models show phase 
transition from free flow to jam in a fundamental di- 
agram. Moreover, the ultradiscrete model includes 
the Fukui-Ishibashi model in a special case. 



1 Introduction 

There are various models of different levels of 
discreteness to analyze the traffic congestion [1]. 
Macroscopic model is defined by a partial differ- 
ential equation based on fluid dynamics and it de- 
scribes a traffic flow by the motion of continuous 
media. For example, Musha and Higuchi used the 
Burgers equation to describe a fluctuation of traffic 
flow [2]. 

System of ordinary differential equation (ODE), 
coupled map lattice (CML) and cellular automa- 
ton (CA) are often used as microscopic model to 
describe each vehicle motion directly. About ODE 
models, time t and vehicle position x are continu- 
ous, and vehicle number k is discrete. CML is sim- 
ilar to ODE but time is discretized [3]. All depen- 
dent and independent variables are discrete for CA 
models. For example, Nagel-Schreckenberg model 
[4], elementary CA of rule number 184 (ECA184) 
[5], Fukui-Ishibashi (FI) model [6] and slow-start 
model [7] are known as effective trafSc model. 
Though evolution rule of CA model is simple due 
to its discreteness, a mechanism of congestion for- 
mation is presented sharply. 




Figure I: Circuit and vehicles. 

Bando et al proposed a noticeable ODE model 
[8]. The model is now called 'optimal velocity 
model' (OV model) and is defined as following. As- 
sume a finite number of vehicles moving on a one- 
way circuit of single lane as shown in Fig. [T] The 
length of the circuit is L and total number of vehi- 
cles is K. Introduce a one-dimensional coordinate 
along the circuit with an appropriate origin. Define 

(t) by a position of vehicle with vehicle number k 
{k = 1, 2, • ■ • , K) at time t. The vehicle number is 
given sequentially to each vehicle as the preceding 
one has a larger number. Note that the preceding 
vehicle k — K is k = 1. Then the evolution 
equation on Xk{t) is 



Xk = A{V{xk+i - Xk) - ife}. 



(1) 
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where A is a constant representing a driver's sensi- 
tivity and V{Ax) is an optimal velocity represent- 
ing a desired velocity of a driver with a distance Ax 
between his vehicle and the vehicle ahead. The ac- 
celeration of fc-th vehicle is determined by ^ and 
is proportional to the difference between its optimal 
velocity and its current real velocity Xk- 
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The typical profile of optimal velocity is shown 
in Fig. [21 This profile reflects a driver's behav- 
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Figure 2: Typical profile of optimal velocity. 

ior; if the distance from the vehicle ahead is short 
(long), he wants to keep low (high) speed. When 
the distance becomes long enough, he wants to keep 
a speed limit of the road. The results obtained by 
the optimal velocity model agree with real traffic 
data well. 

Nishinari and Takahashi reported an interesting 
relation between the Burgers equation and ECA184 
[9]. They proposed a difference equation called 'dis- 
crete Burgers equation' and showed that the Burg- 
ers equation and ECA184 are obtained by contin- 
uum and ultradiscrete limit respectively from the 
discrete Burgers equation. Ultradiscretization is a 
method utilizing a non-analytic limit defined by the 
following formula [10]. 

lim £ log(e^/'' + e^/'' + • • • ) = max(A, B, ■ • • )■ 

(2) 

We obtain a piecewise-linear type of equation called 
'ultradiscrete equation' by ultradiscretizing a differ- 
ence equation. There is a correspondence between 
basic operations of difference equation and of ul- 
tradiscrete one. Usual operations +, x and / of 
difference equation correspond to max, -f and — of 
ultradiscrete one respectively. Thus we can make 
'analytic evaluation' for ultradiscrete equation as 
we do for difference equation. 

Moreover dependent variables can be discretized 
using appropriate initial data and constants of ul- 
tradiscrete equation. Therefore ultradiscrete equa- 
tion is a completely discretized equation in this 
sense. Utilizing this feature, we can show that 
ECA184 originally defined by a binary table is 



equivalent to ultradiscrete Burgers equation. Thus 
asymptotic behavior of solutions to ECA184 can be 
proved by the analytic evaluation reflecting that of 
Burgers equation. As seen by this example, ultra- 
discretization gives a direct relation between CA 
and differential equation via difference one and pro- 
poses a new perspective for CA which can not be 
obtained if we make a closed analysis. 

In this letter, we propose a difference equation 
relevant to the OV model and call it 'discrete OV 
(dOV) equation' . If we take a continuum limit for 
this equation, we obtain (jT]) with a specific V{Ax). 
If we take an ultradiscrete limit, we obtain 'ultra- 
discrete OV (uOV) equation' including ECA184 or 
FI model in a special case. Since uOV equation is 
of second-order on time difference, it can express an 
acceleration effect. Both dOV and uOV equations 
show a phase transition from free flow to jam. 

2 Discrete Optimal Velocity 
Model 

Let us assume the same situation as of OV model 
([1]). The only difference is that a time variable is 
discrete. Assume a time step denoted hy n {n — 0, 
1, •■•) and an interval of time step by 5 (> 0). 
Using these notations, dOV equation is defined by 

- 2x1 + - M log(l + S'V{x-,+, - x^)) 

-log(l + (S(e"^-<"' - 1))}. 

(3) 

If 1 + S'^V{xl^^ - x^) or 1 + Sie'^^-^ ' - 1) in 
the logarithmic terms is or negative, Q is not 
well-defined. However, if 5 is small enough and if 
V{Ax) and initial data are appropriately defined, 
we can easily exclude this problem. 

Replacing a;^ by Xk{nS) and assuming S ^ 0, we 
obtain the following expansion. 

Xk = A{V{xk+i-Xk)-Xk} + -{ik-iif)S+0{5'^). 

(4) 

Thus ([T]) is derived from ([3]) by the continuum limit 
(5-^0 and ([3]) is a discrete analogue to U]). Consid- 
ering this relation, V{Ax) is required to have the 
profile roughly shown in Fig. [H Moreover if we as- 
sume ^ can be ultradiscretized, V{Ax) is required 
to have a more specific form. To realize both con- 
tinuum and ultradiscrete limit, we fix the following 
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form for V^(Ax), 



V{Ax) = a(- TT-r r - - 

where a, b and c are all positive constants. Figure[3] 
shows an example of profile of V{Ax). Figured] 



length and flow Q is equivalent to a total momen- 
tum of vehicles per unit length. Both are defined 
(5) by 




Figure 3: Profile of V{Ax) defined by ^ with a = 
2, 6 = 4, c = 2. 

shows an example of orbits of vehicles. Initial posi- 
tions of vehicles are set at nearly regular intervals 
with small disturbances. Coalescence of jams oc- 
curs at earlier time and three major jams survive 
in this figure. Though not shown in this figure, 
more coalescences occur after a long time. 
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Figure 4: Example of orbits of vehicles for L — 50, 
K = 25,S^ 0.1, A = 1, a = 2, 6 = 4 and c = 2. 

The fundamental diagram is shown in Fig. [5] [1] . 
This diagram shows a dependence of flow Q on den- 
sity p. Density p is a number of vehicles per unit 



p = —(number of vehicles), 
L 



Q 



1 



(ni - no + 1)LS 



m K 



(6) 



n=no k=l 

We can observe three phases, that is, (a) free flow 
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Figure 5: Fundamental diagram with L = 50, 6 = 
0.1, A = 1, a = 2, 5 = 4, c = 2, ?io = 90000 
and Til = 100000. Plotted points are obtained for 
1<K<50. 

phase in a low density region, (b) jam phase in a 
medium density region and (c) tight jam phase in a 
high density region. Since these phases can be ob- 
served for the OV model ([T]), we can consider that 
discretization of time variable in the dOV model 
Q works well. 



Ultradiscrete 
locity Model 



Optimal Ve- 



The dOV equation ^ with optimal velocity Q 
can be ultradiscretized. Let us introduce transfor- 
mation of variable and constants including a new 
parameter e defined by 



+ nS 



-5/e 



(7) 
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Substituting the transformation into ([3]) and ([5]), 
we obtain 



2x1 



„n-l 



= A 



{£log(l + - 



elogfl + e 



(^fc-a^fc )/e 



(8) 

If a, 6, c, (5 are all positive and a < be, we obtain 
the following ultradiscrete equation taking a limit 
e +0. 



^n+l 2r^ri^ 



= yl{max(0, a - max(0, -b{x'^_^_^ - - c))) (9) 
-max(0,4' -4'"^)}. 
Moreover this equation is equivalent to 



2x1 



= A{V{xl+, - xl) - max(0, x^ - xl~^)}, 
where 
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Figure 7: Example of orbits of vehicles for (|T0|) with 
(10) L = bO, K = 2h, A = 0.5, a = 1.9, 5 = 4, c = 3. 



V(Ax) = max(0,6(Aa;-c)+a)-max(0,6(Ax-c)). 

(11) 

Note that 5 is not included in (fTU]) due to the trans- 
formations of xj! and a in ([7]). 

We show a typical profile of V{/S.x) in Fig. El 
Figure [7] shows an example of orbits of vehicles. Po- In this section, we discuss two special cases of uOV 

model. 



4 Special Case of Ultradis- 
crete Optimal Velocity 
Model 




Figure 6: Profile of ^(Aa;) in ([TT|l . 

sitions of vehicles are random integer at initial time 
step. However xj! is generally non-integer since A 
and a are not integer in this example. A fundamen- 
tal diagram using the same constants other than K 
is shown in Fig. [SI Surprisingly three phases clearly 
exist as in Fig.[5l Note that numerical experiments 
are executed by double precision calculation of C 
program. 



4.1 No Overtaking 

Assume A, a, b, c in (fTU)) and (fTTj) are all positive. 
Then if AV{Ax) < Ax, we can derive 



n+i ^ n . n 



Amax(0,a;'fc' 



(a) 

(12) 

Moreover if ^ > 1, (a) is always or negative. 
Therefore we get x^'^^ < x^j^^ on these assump- 
tions. And furthermore, if velocity of all vehicles 
is non-negative, overtaking does not occur. When 
we use the OV (dOV, uOV) model as a numerical 
simulator of concrete traffic flow, overtaking of ve- 
hicle can not occur in a one-way circuit of single 
lane. Though we can avoid overtaking by choosing 
appropriate constants and initial data, assurance of 
no overtaking is important for a real application. 
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Figure 8: Fundamental diagram for (fT(I|) with L — 
100, A = 0.5, a = 1.9, 6 = 4, c = 3, 7io = 1000, 
Til = 2000. Plotted points are obtained for 1 < 
K < 100 and 50 trials are executed for every K. 

4.2 Cellular Automaton 

If constants A, a, 6, c and initial position are all 
integer, any xjj calculated by (fTU|) is also integer. 
Therefore dependent and independent variables in 
((TU)) are all discrete in this case. Moreover, if we 
set A = 1, (fTUl) reduces to 

=2:^+^2:^+1 -4^)+min(0,x^-xr') (13) 
Let us assume V{Ax) > as in Fig. [HI Moreover 



if xl 



> for 



last term min(0, 
= Vixl 



any k 

n-l 



k+1 ' 



at a certain n, the 
in (fT^ becomes and 
> for any k. Therefore 



any vehicle does not go backward if initial velocity 
of any vehicle is not negative. Under this condition, 
(|13p again reduces to the first-order equation. 



V{3 



k+1 



(14) 



Moreover let us consider the case of A = 1, a = 

Wmax, b = 1, C = -ymax + 1 whcrC Wmax IS pOSitivC 

integer. Then V{Ax) in becomes 

V{Ax) = max(0, Ax — 1) — max(0. Ax — Umax — 1) 
'O (Aa;-1<0) 

= < Aa: (0 < Ax - 1< Umax) 
?^max (Wmax < Ax ~ 1) 

(15) 

Assuming a size of vehicles is a unit cell Size, Xy^. I 2^ 
x)! — 1 is a distance between fc-th and (fc + l)-th 
vehicles at time step n. Therefore every vehicle 



moves forward by its distance up to Wmax- This 
model is nothing but the FI model and ECA184 
for Umax — 1. We note that an analogy between 
OV and FI models is commented by Nishinari [11]. 

5 Concluding Remarks 

We propose a new discrete OV model with a dis- 
crete time. Continuum limit of this model is equiv- 
alent to the OV model. We show that orbits of 
vehicles and the fundamental diagram agree with 
those of OV model qualitatively. Moreover this 
model has an ultradiscrete limit and a piecewise- 
linear type of evolution equation is obtained. We 
show that the ultradiscrete OV model also gives 
phase transition in its fundamental diagram by a 
numerical calculation. It includes the FI model as 
a special case. 

We only show a definition, a few features and 
some numerical results about dOV and uOV mod- 
els in this letter. Detailed analysis using various 
combinations of constants is necessary to under- 
stand a dynamics of the models fully. Compari- 
son with other models and with real data is also 
necessary. These points are future problems to be 
solved. 
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